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Abstract 

The  long-range  Interaction  of  two  adsorbed  atoms  mediated  by  the  elastic 
distortion  of  the  substrate  is  calculated  classically  for  an  elastically 
isotropic  substrate.  For  Identical  atoms,  the  interaction  is  repulsive; 
for  different  atoms,  it  can  be  repulsive  or  attractive.  It  varies  as 
0~^  with  the  distance  p between  the  two  adsorbed  atoms.  This  is  the  same 
spatial  dependence  as  for  the  dipole-dipole  interaction  between  two  adsorbed 
atoms.  For  two  xenon  atoms  adsorbed  on  silver,  the  elastic  interaction  is 
about  one  order  of  magnitude  smaller  than  the  dipole-dipole  interaction. 

The  interaction  energy  is  inversely  proportional  to  the  shear  modulus  of  the 
substrate,  so  that  it  may  become  quite  large  near  a distortlve  phase 
transition. 
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I.  Introduction 
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Atoms  adsorbed  on  the  surface  of  a solid  interact  with  each  ocher 

directly  as  well  as  indirectly  through  the  substrate.  The  dipole-dipole 

_3 

interaction  between  two  adatoms  on  a metal  substrate  varies  as  p with 

the  distance  p between  them^^^.  The  oscillatory  interaction  mediated  by  the 

conduction  electrons  in  a metal  was  discussed  by  Grimley,  Einstein  and 

Schrieffer and  appears  to  have  been  experimentally  observed  by  Tsong.^^^ 

We  shall  show  in  a forthcoming  paper  Chat  this  Interaction  decays  more  rapidly 
_3 

than  p . In  the  present  paper  we  Investigate  the  classical  long-range 
interaction  between  two  adatoms  mediated  by  the  elastic  distortion  of  the 
substrate.  Different  aspects  of  the  lattice-mediated  interaction  has  been 
studied  by  Schick  and  Campbell, and  by  Cunningham  et  We 

shall  discuss  these  in  Sec.  VI. 

An  expression  for  the  strain  of  a solid  acted  on  by  external 
forces,  due  to  the  presence  of  an  adsorbed  atom,  is  derived  in  Sec.  II,  and 
an  expression  for  the  interaction  energy  between  two  adatoms  is  obtained  In 
Sec.  III.  A rough  estimate  for  the  magnitude  of  this  interaction  is  also 
given.  Because  of  the  reduced  symaetry  near  a surface,  a quantum  mechanical 
calculation  of  this  interaction  is  rather  complicated  and  was  not  carried  out. 
However,  in  the  appendix,  we  show  that  in  the  case  of  two  Impurities  in  an 
infinite  lattice,  quantum  mechanical  and  classical  calculations  give  the 
same  asymptotic  Interaction  energy,  and  we  expect  that  the  same  is  true 
in  the  vicinity  of  a surface.  We  show,  in  Sec,  TV,  that  the  interaction 
energy  between  identical  adatoms  is  repulsive.  Also  the  interaction  between 
any  two  physisorbed  atoms.  Identical  or  not,  is  shown  to  be  always  repulsive. 
In  Sec.  V,  we  calculate  the  elastic  interaction  energy  between  two  xenon 
atoms  adsorbed  on  the  surface  of  silver,  and  compare  it  with  the  dipole-dipole 
Interaction  of  the  same  system.  A concluding  discussion  is  given  in  Sec.  VI. 
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II  Strain  of  the  Lattice  due  to  an  Adatom 

Consider  an  atom  on  the  surface  of  a solid  occupying  the  space  z^. 

The  adatom  exerts  forces  on  the  substrate  atoms  at  positions  = 

^*i’^i’*i^‘  strain  of  the  lattice  at  a point  ir  far  away  from  the 

adatom,  where  the  direct  interaction  between  the  adatom  and  the  substrate 

atoms  is  negligible,  is  calculated  as  follows:  We  treat  the  substrate  as 

an  elastically  Isotropic  half-space  with  Lame  moduli  g,  and  X,  and  subject 

to  external  forces  acting  at  the  points  The  displacement  u “ (u,v,w) 

at  r = (x,y,z)  in  a half-space  zSQ  acted  on  by  an  external  force  F*(F  ,F  ,F  ) 
~ /g)  - X y z 

at  the  origin  is  given  by 


u 


4nvi  l^r  r3  ^+x  ^r-z 
Attu  [rT  r(r-z)2j 

f is_  r«  _ * — 1 

4ttu  I ^ W.+X  r(r-z)  I 
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r(r-z 
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(2.1) 


with  a similar  expression  for  v.  The  z-component  of  tbs  displacement  is 
not  needed  here.  We  now  sum  the  displacements  due  to  the  forces 
exerted  by  the  adatom  at  the  positions  The  forces  are  assumed  to  be 

central.  For  a system  with  3-fold  or  4-fold  symmetry  about  the  axis  through 
the  adsorbate  normal  to  the  substrate  surface,  and  taking  into  account  that 
the  forces  on  the  substrate  atoms  must  add  up  to  zero,  the  displacement  at 
a point  near  the  surface  far  from  the  adatom  is  found  to  be 


“ “ ^ “r^  1 

I (2.2) 


r 


- "V.  — 
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with  = (j^  ,Yi  ,0)  • 

III  Interaction  Energy  of  Two  Adatoms 

The  energy  of  the  system  consisting  of  the  lattice  and  an  adatom  relative 
to  that  of  the  unstrained  lattice  is 


E = - 

i 


(3.1) 


For  two  adatoms,  (1)  and  (2)  at  and  £2  respectively,  (see  Fig  1),  the  ' 
elastic  energy  is 


i 

where  is  the  force  on  the  substrate  atom  i exerted  by  adatom  a 

and  u|**^  is  the  displacement  of  atom  i due  to  F? 

The  interaction  energy  is  given  by 


(3.2) 


int 


,(» 


(3.3) 


where  ^12  “ ^2  " £.1 

mrwi  c “ - is  the  Poisson  ratio  of  the  substrate. 

2 (n+X) 


The  interaction  energy  can  also  be  expressed  in  terms  of  the  strains 
produced  in  the  subatrate  by  the  adatoms.  We  introduce  the  two-dimensional 
dilatation  (due  to  adatom  o ) 


(3.4) 
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The  quantity  - o p is  ^ constant,  independent  of  position.  Substitu- 

ting this  into  (3.3)  we  obtain 

^Int  ■ 

To  obtain  an  order  of  magnitude  estimate  for  the  strength  of  the 
interaction,  let  us  assume  that  at  |p-p  [=5a.u.  (which  corresponds  roughly 

Of  ^ k 

to  the  position  of  the  substrate  atom  nearest  to  the  adatom)  0^  0.1. 

_3 

For  most  metals  O sa  h snd  ^ 10  a.u.  This  gives,  for  two  adatoms 

10  a.u.  apart,  an  interaction  energy  of  about  O.OSeV. 

IV  The  Sign  of  the  Interaction  Einergy 

The  interaction  between  two  adatoms  may  be  repulsive  or  attractive, 
depending  on  the  relative  signs  of  the  sums  in  parentheses  in  (3.3), 
or  of  0^^^  and  in  (3.4). 

For  identical  adatoms,  they  have  the  same  sign  and  hence  the  Inter- 
action is  repulsive. 

We  shall  now  show  that  It  Is  also  repulsive  between  physisorbed 
atoms,  where  the  repulsive  potential  between  the  adatom  and  the  substrate 
atoms  is  of  shorter  range  than  the  attractive  potential.  Since  the  adatom 
is  in  equilibrium  we  have,  for  forces  in  the  z-directlon 

r (?<■•>- ^-0  ».i) 

where  end  F^  ^ ere  the  megnleude*  of  the  repulsive  end  attrecclve 

parts  of  the  force  IF^  respectively  and  Ri  Is  the  distance  between  the 
adatom  and  the  substrate  atom  at  site  1,  with  the  z-component 
(see  Fig  2). 

(r)  (a) 

Since  the  repulsive  force  is  of  shorter  range,  we  have  F 2 F for 

n n 

sites  n within  a certain  distance  from  the  adatom  and 

F(r)<  p(i)  £qj.  B farther  away 

mm 
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(4.1)  can  be  written  as 


E f 


E <^1”’  - r 

m m 


(4.2) 


with  each  term  on  either  side  of  the  equality  positive. 


Now,  for  each  layer  of  substrate  atoms,  with  constant  z 


and 


Y'  ^p(a)  _ p(r).  ^ ^ y /p(a)  _ 


m 


so  that 


E (f'">  - fi‘>) 


(4.3) 


We  now  examine  the  sum  in  (3.3), 


E(Fl-£j)-E<-{'’-fJ*’)R7 


E ^ r E ffi*’  - 


i, 

?SL 

R 

m 


«#W  a ‘A, 
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Ro 


ErF'o' - < 0 

■'  " y mm  my_ 


m 


m 


(4.4) 


Thus  (4.4)  is  always  negative,  l.e.,  the  Interaction  energy  (3.3)  is 
repulsive  between  physisorbed  atoms. 


V Application  to  Xenon  Atoms  Adsorbed  on  Silver  Surface 

As  a specific  example,  we  consider  the  adsorption  of  noble  gas  atoms 

on  the  (111)  face  of  a fee  crystal,  with  the  adatoms  at  centered  sites 

(synmetrlcal  between  three  surface  atoms).  The  attractive  potential  is 

taken  to  be  of  the  fonn-jyJ^  ^ and  the  repulsive  part  taken  to  be 

i 

some  short-range  potential.  We  make  the  approximation  that  silver  is  an 
isotropic  material. 

If  we  consider  only  the  attractive  interaction  between  the  adatom  and 
the  three  nearest  substrate  atoms,  we  have 


where  a = distance  between  two  nearest  neighbor  substrate  atoms 

Rq*  distance  between  the  adatom  and  one  of  the  nearest  substrate 
atoms 

Fq"  component  parallel  to  the  substrate  surface  of  the  attractive 
force  between  the  adatom  and  one  of  the  nearest  substrate 
atoms. 


Including  the  repulsive  interaction  and  the  interaction  between  the 
adatom  and  other  substrate  atoms  in  the  sum,  we  obtain  the  value 

^ * -^0  “ * (5-1) 


J 


1 
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where  cy  Is  a dimensionless  number  which  depends  on  the  distance  of  the 
adatom  from  the  substrate  surface. 

The  interaction  energy  is  thus 


int 


9a" 

2nti 


(1-  a) 
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Aa 


4 1-a 
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■y  -3 


2ttm,W(Z-Z  )V  P 
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(5.2) 


where  the  attractive  potential  is  now  written  in  the  form^^^^(see  Fig  3) 


(Z-Zo)3 


(5.3) 


with  Z the  position  of  the  adatom  relative  to  the  first  layer  of 
substrate  atoms 

Zq  the  position  of  the  "reference  plane" 

C the  van  der  Waals  coefficient 
and  A a dimensionless  constant  depending  on  a 


For  xenon  adsorbed  on  silver,  we  have,  in  atomic  units,  C > 0.813, 


Z - 6.85,  Z-Zq 


4.10<“), 


a “ 5.47, 
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0.6  and  A ■ 0.65. 


This  gives,  for  xetion  atoms  separated  by  a distance  of  10  a.u.  (which 

corresponds  to  monolayer  coverage),  an  Interaction  energy  of  roughly 

-4  -3 

10  eV.  The  interaction  is  repulsive  and  varies  as  p with  the 

distance  P between  the  two  adsorbed  atoms.  This  is  the  same  spatial 

dependence  as  for  the  dipole>dipole  interaction  of  the  two  adatoms,  which 


is 


(1) 


E 


dlp<^ip 


-3 
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where  and  V-2  are  the  dipole  moments  of  the  adatoms  1 and  2 . 

For  xenon  atoms  adsorbed  on  silver  at  monolayer  coverage,  p,  = 0.2  debye 
giving 


(11) 


lO’^eV. 


VI  Conclusion 

The  elastic  interaction  energy  of  two  adatoms  is  repulsive  for  identical 
atoms  and  also  for  any  two  physisorbed  atoms.  For  arbitrary  adatoms,  the 
interaction  may  be  attractive  or  repulsive.  The  spatial  dependence  is 
p“^,  which  is  the  same  as  for  the  dipole-dipole  interaction.  For  xenon 
adsorbed  on  silver,  the  elastic  interaction  is  roughly  one  order  of  magnitude 
smaller  than  the  dipole-dipole  Interaction  energy,  and  the  same  is  true 
for  xenon  on  other  noble  metals.  For  chemisorbed  atoms,  the  interaction 
energy  between  pairs  of  adatoms  at  monolayer  coverage  is  roughly  0.01  to 
0.1  eV.  We  note,  from  (3.3)  that  the  interaction  energy  is  inversely 
proportional  to  the  shear  modulus  p,  of  the  substrate,  so  that  it  may 
become  quite  large  near  a distortive  phase  transition. 

Finally,  we  would  like  to  compare  our  work  with  two  previous  calcula- 
tions of  the  interaction  energy  between  adsorbed  atoms  mediated  by  the 
substrate  lattice. 

Schick  and  Campbell  calculated  the  phonon-mediated  interaction  energy 
of  non-localized  helium  atoms  adsorbed  on  argon-plated  copper.  They 
treated  the  adatoms  as  Bloch  waves,  with  two-dimensional  quasimomentum 
parallel  to  the  surface,  in  contrast- to  the  present  calculation,  wViich 
deals  with  atoms  localized  at  definite  sites.  Cunningham  et  al.  calculated 
the  total  phonon  zero  point  energy  of  a simple  cubic  lattice  with  two 
adatoms  localized  at  two  sites  a distance  p apart.  They  obtained  a 
weak  attractive  interaction  behaving  asp  at  large  distances.  This 
interaction  may  be  regarded  as  the  leading  quantum  correction  to  the 
classicaldlstortive  interaction  calculated  in  the  present  paper  * (In  the 
particular  oversimplified  model  of  Cunningham  et  al.  the  classical  distortion 
is  exactly  zero).  For  adatoms  at  nearest  neighbor  sites  with  all  masses 
and  direct  interactions  taken  as  comparable  they  find  for  their  indirect 
interaction 


E , - - 10 

zero-point 


hu). 
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where  is  a characteristic  frequency  of  the  lattice.  With  the  same 
assumptions  the  classical  distortion  energy  calculated  in  the  present  paper 

-2  2 

distortion  o 

where  K is  an  interatomic  force  constant  and  a is  the  lattice  constant. 

o 

Typically,  then, 

-3  -4 

E /E,.  . ~ 10  or  10  (nearest-neighbors) 

zero-point  distortion 

t'.ecause  of  the  more  rapid  fall-off  of  the  zero-point  interaction  compared 
to  the  distortive  interaction,  the  ratio  becomes  even  smaller,  for  more 
widely  separated  adatoms. 

One  of  us  (K.H.L.)  acknowledges  with  thanks  the  award  of  an  IBM 
graduate  fellowship. 
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APPENDIX:  Elastic  Interaction  Energy  of  Two  Defects  in  an  Infinite  Lattice 

We  calculate,  both  quantum  mechanically  and  classically,  the  elastic 
energy  of  a large  isotropic  crystal  of  volume  n with  two  defects  in  the 
interior,  at  r.and  r_  , which  exert  forces  (per  unit  volume)  ^ a.  i 


f(2) 


on  the  atom  at  sitd  /and  shall  verify  that  the  results  are  identical 


for  large  separations. 

To  calculate  the  energy  quantum  mechanically,  we  write  the  perturbing 
potential  as 


and  express  the  displacement  vector  in  terms  of  phonon  operators 


v|tj(  Oj  4<-a0 ~ £j  <1> 


with  wavevector  £,  polarization  index  j and  polarization  vector 
We  express  the  force  in  terms  of  the  polarization  vectors 

= f(1)  + f^2)  {Fj^^(q)e^3'^-"-l^  0^(5)  + ~ -2^(3)} 


a.j 


The  energy  of  the  system  relative  to  its  ground  state  is,  by  second-order 
perturbation  theory. 


E 


|(0l6Vlq,j)l^ 

-hoJj 


(q)Fj^^  C-q) 

a'j  (q^) 


cos 


Here  d is  the  density  of  the  crystal,  and  R « rj"  The  second  term 

on  the  right  side  of  (A. 4)  represents  the  elastic  interaction  energy  of 


(A.l) 


(A.2) 


(A.  3) 


J 
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the  defects; 


•"int  - “ d jf  j a^'TqT~~  <3-  ' 3. 


(A.S) 


The  classical  elastic  energy  of  the  system  is  given  by 


- - Tf  (r)  ■ u(r)d^r 


(A. 6) 


with  the  displacement  u(r)  determined  by  the  equation 


u7  u + (u.  + X)  7(7-u)  + F = 0 


(A.  7) 


Taking  the  force  F(r)to  be  of  the  form  (A. 3)  but  with  the  label  t replaced 
by  the  continuous  variable  r,  the  solution  for  (A, 7)  is  found  to  be 


u(r)  = 


u 

A 

V 

\j 


f“  fV  f“ 

/ Ak-js inflcos f,7  + 1 cosScoscs  — *^2  sin 

^2p+\  ^ H ' IT 


=Ea 

a <1 


F?  - - . F^,  _ . Ft2 


^ sin0sinc?-j — - cosQsincp  + — — coscp 
2n+X  M.  H 


,iq* (£-£«) 


' — ii—  COS0  — _£l  sing 

2p,"^X  p, 


/ 


/ 


(A.  8) 


with  0 and  cp  the  polar  and  azimuthal  angles  of  q 
a = 1 or  2 denoting  the  defects  I and  2 
and  j = L,  T]^  or  T2  referring  to  the  longitudinal  or  transverse  components. 

The  elastic  energy  is  thus  E =•  - £, 


2u+X 


■)- 
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oE- 

q 


cos (q -R) 


■) 


The  interaction  portion  o£  this  energy  can  be  rewritten  as 


cl 

E . 
int 


2 2 
Cj(0)q^ 


cos  (q*R)d^g^ 


(A.  9) 


(A. 10) 


where  ^0)  the  (isotropic)  sound  velocity  at  q = 0 . 

We  note  that  (A. 10)  is  the  same  as  (A. 5)  except  for  the  presence  of 

2 2 1 2 
Cj(0)q  in  the  denominator  in  place  of  “>j  (q^)  . ^ shall  now  show 

that,  for  large  R,  the  main  contribution  to  the  integral  in  (A. 5)  comes  from 

qadO,  so  that,  asymptotically,  the  two  expressions  for  the  elastic  interaction 

energy  are  identical. 

The  integrand  in  (^.5)  is  a function  only  of  the  magnitude  of  q and  the 
angle  S between  j|  and  R .It  can  be  expanded  in  a Fourier  series  in  S,  or  in 
a power  series  of  the  form 


E ^ (q)F (-q)cos  (q*R)=  E(a  (q)cos"9'H>  (q)s in9cos*^9) 
- „n  n 


(A. 11) 


so  that  (A, 5)  becomes 
2 

E,  = - — j 1 E(a  (q)cos'^0  +b  (q)sin9cos"0)  d:os9 

int  <j  J J n n 


(iR)" 


(iR)" 


For  large  R,  we  find 


(12) 


J ^n 


O 


dq 


(A.12) 
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.n 


,2k-l  , , 

. - , n I i'-^  ^ -I  , , iqR  , 1^)1  **  3"  ^ sin  qR  'Nj 

.(.1^  K^m  , —IT— ‘I  k L ^ 

-00  dq  dq  ()q  ^ q=0 


d"a 

= (.1)"  !^ 1. 

^ 2R  . n 
dq 


.2k-l 

> ^ 0 3 (q),n-2k  . . „ ^ 

q=0  R k dq  dq 


qR  yj 


q=0 


(A. 13) 


where  H(q)  is  the  Heaviside  function 


H(q)  = { 


0 q < 0 


(A. 14) 


1 q ^ 0 

Similarly,  for  m = n or  n + 2 (cf.Eq.  (A. 12)) 


O oq 


p d"h)  (q)  nr-  (q)  ,m-2k 

<-')"  j + ^[■V7w:r'  -ST2T  V’MJ  „ 

O dq  ’ ^ k dq  oq  q=0 


d^*^"V(q)  3®-2k 


UH 


a=0 


./rd% 


dq“ 


.2k-l, 


^ ‘^<-2'  [— 2CT-  ^2k  . 

q-0  R dq  dq 


q*0 


(A, 15) 


Thus  the  main  contribution  to  the  integral  in  (A. 5)  comes  from  «f!dO,  which 
is  what  we  wanted  to  {.rove. 

This  establishes  the  Identity  of  the  classical  and  quantum  mechanical  ' 
result . 
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ho  lonK-range  interaction  of  two  adsorbed  atoms  mbdlated  by  the  elastic  distortion 
of  the  substrate  is  calculated  classically  fo^.  an  elastically  isotropic  substrate. 
For  identical  atoms,  the  Interaction  is  rejjdisive;  for  different  atoms,  it  can  be 
repulsive  or  attractive.  It  varies  as  p”3  with  the  distance  p^b^tween  the  two 
adsorbed  atoms.  This  is  the  same  spatial  dependence  as  for  the  dipole-dipole  inter- 
action between  two  adsorbed  atoms.  For  two  xenon  atoms  adsorbed  on  silver,  the 
elastic  interaction  is  about  one  order  of  magnitude  smaller  than  the  dipole-dipole 
interaction.  The  interaction  energy  is  inversely  proportional  to  the  shear  modulus 
of  the  substrate,  so  that  it  may  become  quite  large  near  a distortive  phase 
t ransition . 
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